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This paper discusses various wormhole solutions in a spacetime with an equation of state that is 
both space and time dependent. The solutions obtained are exact and generalize earlier results on 
wormholes supported by phantom energy. 
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INTRODUCTION 



Traversable wormholes, whose possible existence was 
first conjectured by Morris and Thorne in 1988 J\, are 
actually shortcuts that could in principle be used for 
traveling to remote parts of our Universe or to different 
universes altogether. The meticulous analysis in [l| has 
shown that such wormholes can only be held open by the 
use of exotic matter. Such matter violates the weak en- 
ergy condition (WEC), which requires the stress-energy 
tensor to obey Tap^i^fi^ > for all time-like vectors 
and, by continuity, all null vectors. For example, given 
the radial outgoing null vector (1, 1, 0, 0), we obtain 



T, 



= P + P>o. 



(Recall that T^^ — p and Tff = p in the orthonormal 
frame of reference.) So if the WEC is violated, we have 
p + p<0. 

Interest in traversable wormholes has increased in re- 
cent years due to an unexpected connection, the discov- 
ery that our Universe is undergoing an accelerated ex- 
pansion 0, S]- This acceleration, caused by a negative 
pressure dark energy, implies that a > in the Fried- 
mann equation ^ ~ — ^(p + ip). (Our units are taken 
to be those in which G = c = 1.) The equation of state is 
p = —Kp, if > i, and p > 0. While the condition K > ^ 
is required for an accelerated expansion, larger values for 
K are also of interest. For example, K = 1 corresponds 
to a cosmological constant. 

Of particular importance for us is the case K > 1, 
referred to as phantom energy. For this case we have 
p -\- p < 0, in violation of the weak energy condition. 
As noted realier, this condition is the primary prerequi- 
site for the existence of traversable wormholes. (Strictly 
speaking, the notion of dark or phantom energy applies 
only to a homogeneous distribution of matter in the Uni- 
verse, while wormhole spacetimes are necessarily inho- 
mogeneous. However, the extension to spherically sym- 
metric inhomogeneous spacetimes has been carried out 

In a recent paper, F. Rahaman et al. [5| discussed 
wormhole solutions that depend on a variable equation 



of state, i. e., ^ — —m{r), where m(r) > 1 for all r. 
The variable r refers to the radial coordinate in the line 
element 



.2/(r) 



1 



-dr^ +r'^{d9'^ + sin'^ed(j>^); (1) 



in other words, m = m{r) is independent of direction. It 
is shown that given a specific shape function b — b{r), 
it is possible to determine m = m(r) and vice versa. It 
is also assumed that f'{r) = 0, referred to as the "zero 
tidal-force solution" in Ref. 

An earlier study ^ assumed that the equation of state 
is a function of time. In this paper we will assume that 
m — m(r, t) is a continuous function of both r and t. As 
in Ref. 5], however, we retain the assumption that the 
function values are independent of direction. For reasons 
that will become apparent later, we also assume that the 
change in t is very gradual. 

The main goal in this paper is to show that the time- 
dependent metric describes a slowly evolving wormhole 
structure without assigning specific functions to b and m. 
Moreover, the function / in line element llj need not be 
a constant. 

Evolving wormhole geometries are also discussed in 
Refs. and 0. 



II. THE PROBLEM 

Normally, one would begin with the general line ele- 
ment 

ds' = -e^^Wdt^ +g2aM^^2 ^^2(^^2 ^gjj^2^^^2)_ (2) 

In view of line clement ([T]) , 

2a(r) _ 



1 



1 — b{r)/r 

As already noted, b — b{r) is the shape function; fe(ro) — 
To, where r = rg is the radius of the throat. As a result, 

lim a{r) = +oo. 

r—>ro + 

Rccall that "f{r) is referred to as the redshift function. 
This function must be finite everywhere to avoid an event 
horizon. 



2 



Returning to the function m = m(r,t), since, for any 
fixed the function is invariant under rotation, we need 
a time-dependent metric with the same property 

ds^ = ~e^^^^'''>dt^ + e^"^-'*'>dr^+r^{de^+sm^edcP^). (3) 

This metric describes an evolving wormhole. 

Observe that the shape function is now given by 



(4) 



It will be seen in the next section that a{r,t) depends 
directly on m{r,t). So if m{r,t) changes slowly enough 
with respect to time, then the last term on the left-hand 
side of inequality (O becomes negligible, ensuring that 
the WEC will always be violated. 

From the Einstein field equations G^^ = ^""-^a/? ^'^^ 
the equation of state p — —m{r,t)p, we have G^f = Snp 
and Gff — 8TT[—m{r,t)]p. Using Eqs. (O and ([6]), we 
obtain the following system of equations: 



To obtain a traversable wormhole, the shape function 
must obey the usual flare-out conditions at the throat, 
modified to accommodate the time dependence: 

d 

b{rQ,t) = ro and -g-b{rQ,t) < 1 for all t. 

Another requirement is asymptotic flatness: 6(r, t)/r — > 
as r —^ oo. 

The components of the Einstein tensor in the orthonor- 
mal frame are available from Ref. Q: 



G; 



-2a{r,t) 



d_ 

dr 



1 



a(r,i) + 4,(l-e-W)), (5) 



or 



r 



G,, = 8.T,, = ^e--(-*)|:a(r,t) 



-2a{r,t) 



Gff = SirTff — 87r[— m(r, 



or r^ 



r 



After substituting and rearranging terms, we have 



d 



mir, t)—a{r, t) 



d_ 

dr 



7(r,t)-l(e2"('-'*)-l) [m(r,t)-l]. (10) 



G,-, = ^e-^('-^*)e-"('-'*)^«(r,i), (7) 



Goo — Gil 



5' . ^ 



d d ( d 

--^[r^t)-air^t)+\-a{r^t) 



-1a{r,t) 



Q^2 



j{r, t) 



d d f d 



1 



-2a{r,t) 



d , . d , . 



(8) 



Recall that from the Einstein field equations in the 
orthonormal frame, G^^ = SttT^^, the components of the 
Einstein tensor are proportional to the components of the 
stress-energy tensor. In particular, T^- = T-j- = g^Gj- = 
±/ is interpreted as the energy flux in the outward radial 
direction The WEC now becomes p + p ± 2/ > 0. 
So if the WEC is violated, then 



±^e-'^(''^*)e-"('''*)^a(r, t) 



d_ 



<0. (9) 



III. THE REDSHIFT FUNCTION 



In this section we solve Eq. (|10p by letting the redshift 
function take on a speciflc form. One obvious choice is 
7(r, t) = constant, so that ■^j{r,t) = 0; the other is 

-§-j{r,t) = Eii^iliLL^ allowing the solution of Eq. (jTU]) by 



separation of variables. 

If #7(r, t) = 0, then Eq. ^ becomes 



d 



dr 

Rewritten as 



2|:a(r,t) 1 



g2a(r,t) _i J, y m{r,t) J ' 

one recognizes the form / -r^j ~ In (e" — 1) — u. The 
result is 



In (e2"('^^*) - 1) -2a(r,t) 
Next, we solve for e^"'''^*^: 



-Inr 



dr' 



r'm{r' ,t) 



„2Q(r,t) 



rr dr' 



1 



Evidently, 



b{r, t) = e^^ '■'"•(''■'•*) . 



(11) 



(12) 
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At r = ro, we have 



rr-o dr 
gJc rn.(r,t) 



1, 



(13) 



which imphes that 5(ro, t) = tq for aU t. Also, 

1 



rm{r, t) 
1 



1 



< 1. 



So the flare-out conditions are met without any addi- 
tional assumptions on either b or m. Since m = m{r,t) 
is continuous, the constant c can be uniquely determined 
from Eq. (fT3|) : 



dr 



rm{r^ t) 



In To 



for any fixed t. (So c is actually a function oi t.) 
The line element is now seen to be 



ds^ = -de 



1 - 



dr^ 



+ r'^{de^ +sm^0d(l>^). (14) 
As an illustration, in the special case m(r, t) = K, we 



obtain c ^ 



The result is 



„2a(r) _ 



which is Lobo's solution [ll| . 

Remark: Returning to inequality ([9]), since a(r, i) 
depends on m{r,t), {d/dt)a{r,t) is small only if 
{d / dt)m{r,t) is small, explaining our earlier requirement 
that m(r, t) change only gradually with respect to time. 

For the other choice of 7(r, i). 



d_ 

dr 

Eq. (fTU]) becomes 
or 

m(r, i) — 1 



7(r,t) 



m(r, t) — 1 
Yr ' 



1. |^g2a(r,t) _ [m(r,<)-l] 
^[m(r,<)-l] (-1 - e2"('-'*) + l) , 



2r 2r 
1 

~ 2^ 



2^a(r,t) _ 1/ ^ 

g2Q(r,t) ^ ^ m(r, i) 



It follows that 

g-2a(r,t) 

So by Eq. 61), 

6(r, t) = r ^1 — e 



1 



m(r', i) 



dr'. 



-2a(r,t) 



1 



1 



7n{r' , t) 



dr' 



It now becomes apparent that c — ro for all i, since we 
must have 6(ro , i) = rg . Thus 



5(r, t) = r 



1 



m{r' , t) 



dr' 



(15) 



while the line element becomes 



ds' 



dr^ 



f Ml 

Jrn r' \ 



i(r' ,t) 
q2 I ,-„2 



dr' 



+ r\de^ + &\TL'ed4)^). (16) 



Once again, 
d 



dr ^ ^ rn{ro, t) 



< 1 for all t. 



Returning to the redshift function, we have up to this 
point 



7(r, t) 



m{r\t)_]_ 
2r' 



(17) 



The constant ci will be obtained from the jmiction con- 
ditions, described below. 

Based on previous studies involvin g st atic wormholes 
supported by phantom energy [ll|, fl^ Il3l| , our spacetime 
is not likely to be asymptotically flat. The wormhole ma- 
terial will therefore have to be cut off at some r = a and 
joined to an external Schwarzschild spacetime. Moreover, 
in Ref. [13], b = b{r) actually attains a maximum value 
at some r = a, which then becomes a natural place at 
which to perform the junction. Accordingly, we will as- 
sume that for any fixed t, &(r, t) has a maximum value 
at some r = a. (While other values could be chosen, the 
choice suggested here yields a particularly elegant solu- 
tion.) So we proceed by determining the critical value: 



d_ 
dr 



b{r, t) = 1 



1 



1 



dr' 



m{r' , t) 
r \ m(r, t) 



= 0. 



By assumption, equality holds for some r — a. As a 
consequence. 



a ^ 



1 



1 



m(r, t) 



dr 



1 



m{a, t) 



(18) 
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As noted in Ref. to match our interior solution 

to the exterior Schwarzschild solution 



If 



1--) dr^ 



+ r'{d9' + sin' 9 del}') 



at r — a (t fixed) requires continuity of the metric. Be- 
cause of the assumption of spherical symmetry, jthe com- 
ponents ggg and g^^ are already continuous As a 
result, the continuity requirement has to be imposed only 
on the remaining components. For every t 

ff£t(int) («) = .9tt(cxt) («) and gff (int) (a) = 5ff (cxt) (a) 

for the interior and exterior components, respectively. 
These conditions now imply that (for every t) 



7int(a) = 7oxt(a) and foi„t(a) = &cxt(a)- 



Hence 



^2a{a,t) 



-. _ b{a,t) 1 _ 2M 
r „ n 



We now see that the total mass of the wormhole for r < a 
is given by M = ^b{a, t) for any fixed t. So by Eqs. 
and CH]) 



M ^-a(l 

2 V m{a,t) 



By Eq. 1^, 



1 



2 1 , 

= 1 -all , , 

a 2 V m(a,t) 



or 



=27(a,t) 



m(a, 



(19) 



Eq. ([TO]) can now be used to determine ci = ci (t) , thereby 
completing the line element, Eq. pBl) . 

As an illustration, if m(r, t) = i^T, then a = roe-'-^'^^^^'' 
[from Eq. (HH)] and Ci = ro(ii:e)i/(-^-i). 



IV. ADDITIONAL SOLUTIONS 

As noted in Ref. [13], to obtain additional exact solu- 
tions, 7 must depend directly on a. The corresponding 
condition for the time-dependent case is 

^7(r,i) = i^[a(r,i)]^a(r,i) 

for some elementary function F. Since these cases are 
just extensions of the cases discussed in Ref. [1^, we will 
merely summarize the results. 



F[a(r,t)] - 



m(r, t) 

' p2a(r,t) ' 



then Eq. becomes 
d 

mir, t)—a{r, t) 



= ' } - TT (e^"^'^'*) - 1) [m(r,t) - 1 . 

The solution is 

g2a(r,t) ^ I 

r J, f 1 _ _i ) dr' ■ 

Jro r \ m{r .1) I 

(If m(r, t) = K, the solution reduces to 

2Q(r) _ _ 



1 



discussed in Ref. UM-) Also, 



l-l/K ' 



6(r, t) ~ r 



1 



'ro V m(r',i) 
As in the previous cases, 

d 1 

5(ro,i) = ro and — 6(ro,t) = — — < 1 

or m{ro,t) 

for all t. 

The determination of the redshift function and the 
junction to an exterior Schwarzschild solution follows 
along the lines discussed in Sec. IIIII 

Another solution comes from 

F\a(r,t)] ^ ^ , \,' ' . 

Substituting in Eq. pop and simplifying, we get 

+ 1 ("l - = ^^"("'^) 

e2a(r-,t)_i 2r \ m{r,t)) e4"('-.*)-l' 

The solution is 

-1 



„2a(r,t) 



pJro r' V m(r',t) ) _ \ 



(If m(r, t) = K, this reduces to 

p2a(r) _ _ 



1 



also discussed in Ref. [l^.) Here 

6(r,i) =r(l-e-2"('-.*)^ 

Once again, 

d 1 

6(ro,i) = ro and -^6(ro,t) = ^ -r < 1 

ar m(ro,t) 

for all 
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V. SUMMARY 

This paper discusses several exact solutions of the Ein- 
stein field equations describing traversable wormholes 



supported by a generalized form of phantom energy: the 
equation of state is given hy p/ p = —m(r, t), m(r, t) > 1. 
The function m — m{r,t) is a continuous function of r 
and t, invariant under rotation, slowly evolving in time. 
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